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Abstract 

We complexify one of the Natanzon's exactly solvable potentials in 
VT symmetric manner and discover that it supports the pairs of bound 
states with the same number of nodal zeros. This could indicate that 
the Sturm Liouville oscillation theorem does not admit an immediate 
generalization. 
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I. FRAMEWORK: VT SYMMETRIC QUANTUM MECHANICS 



Bound states in a smooth real potential V(x) > are most easily interpreted 
in the language of Sturm-Liouville theory Its oscillation theorems imply that 
the N— th bound state 4>n{x) possesses N nodal zeros. Even the standard boundary 
conditions may be understood as the presence of an additional pair of zeros which are 
located at both ends of the interval of coordinates J = (—00, 00) in one dimension 
or of the radial axis J = (0, 00) in three dimensions. 

A non-standard, additional constraint is only necessary for some strongly singu- 
lar forces @. Unfortunately, the latter, mathematically rigorous requirement of the 
elimination of the irregular solutions may prove fairly counterintuitive and, for this 
reason, it is often being forgotten in practice. For illustration, we may recollect the 
paper on supersymmetry || which describes a singular map of the shifted harmonic 
oscillator. This example (attributed to A. Khare in Acknowledgements) is marred 
by a subtle violation of the boundary condition in the origin. A very similar incon- 
sistency appears in the "conditionally exact" model by A. de Souza Dutra f|, in 
some phenomenological studies in quantum chemistry || etc. 

An extensive clarification of the latter point may be found elsewhere M. In 
order to avoid similar misunderstandings, the underlying mathematics can be sig- 
nificantly simplified via and analytic continuation (or even re-formulation) of the 
physical boundary conditions J7|. This is an innovative idea which proved unexpect- 
edly fruitful. Recently, one of its more specific versions became a cornerstone of the 
so called VT symmetric approach to quantum mechanics ||. 

The latter formalism dispenses with the Hermiticity of the Hamiltonian H and 
keeps only a weaker requirement of the commutativity of H with the product of parity 
(Vx = —x) and complex conjugation (Ti = —i which mimics the time reversal). The 
singular Schrodinger equations may be then regularized in the way which preserves 
the reality of the spectrum. Such a VT symmetric non-Hermitization can be applied 
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easily to many solvable potentials, e.g., via a constant complex shift of the coordinate 
axis |§, 



r = r(t)=t — ie, t 6 (—00,00), e > 0. (1.1) 
This enables us, in essence, 

• to work, more easily, with the analytic wave functions ip{r) which prove avail- 
able within the whole complex plane of r, 

• to return, whenever necessary, to the Hermitian Hamiltonian by means of the 
limiting transition e — > and/or of a suitable selection of the additional con- 
straints. 

Both the one-dimensional and three-dimensional exactly solvable harmonic oscil- 
lators V {E \r) = £(£ + l)/r 2 + r 2 with parity £ = -1,0 or angular momentum 
£ = 0, 1, . . ., respectively, become suddenly tractable on equal footing (details may 
be found in ref. |1U|| ). In the three-body context of the so called Calogero's exactly 
solvable model |ll], the same singular harmonic-oscillator- like differential equation 
and the same limiting transition e — > reproduce the correct Hermitian solutions 
in spite of the utterly different physical origin and meaning of the singular term 



12]. Via a modified transition to the Hermitian limit, another solvable model is 



revealed [pf] . 



One can summarize that the VT symmetric quantum mechanics offers a new 
approach to the explicit construction of bound states. It will be applied here to the 
potential 

V(f) 3/4 . 2»/?exp(2»fl C 

lU (1-e 2 ^) 2 Vl - e 2i « l-e 2i « { ' } 



with certain unusual features. Already Singh and Devi |14| have noticed that it does 



not belong to the family of the shape invariant potentials. Dutt et al emphasized 
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that this is the only existing nontrivial example of the conditionally exactly solvable 
potential. We have recently shown |16| that this exactly solvable potential of the 



Natanzon class |17| possesses a natural supersymmetric interpretation and offers a 
nontrivial opportunity of the rigorous study of its solutions near the singularity. The 
latter point proved in fact also the main motivation of our continuing interest in this 
example. 



II. THE METHOD: A CHANGE OF COORDINATES 



In a way which dates back to the work of Liouville |y| one can change variables 
in the Schrodinger equation 



(2.1) 



using the recipe 

r - £ = £(r), ^ E \r) -> ^ (D) (0 = ^'[r(0] • ^KO]- 

This generates the new differential equation of the similar form, 



[2.2) 



^- 2 + y (D) (o 



(2.3) 



The explicit relationship between the two respective potentials is determined by the 
identity which is easily derived, 



VW(r)-EW = [e(r)f{v( D \t)-EW} + °~ 
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[C'(r)J 



(2.4) 



Once you postulate the solvability of the original equation (|2.1|) in terms of some 
known polynomials, it is possible to construct another exactly solvable equation by 
the suitable choice of the re-parametrization £(r). Its special cases which mediate 
all the mutual canonical transformations between the real shape invariant potentials 
were listed in review Jl9| (cf., in particular, Figure 5.1 there). 
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Once we switch our attention to the non- Hermit ian examples and retain the split 
of the shape invariant models into the Laguerre-solvable and Jacobi-solvable subsets 



(cf. Figure 5.1 of ref. |19j once more), we can easily parallel many of the above 
transformations. For the Laguerre-solvable subset, all the details may be found in 
refs. and ||21|| . In a less exhaustive manner, the Jacobi-solvable shape invariant 



VT symmetric family has been described in refs. || and . One of the most striking 
features of some of the complexified changes of variables lies in the characteristic 
bent shape of the integration contours £. For illustration, let us recollect the implicit 
definition 



sinhr 



-le 



[2.5) 



of the Hulthen-force generating function £ = £(r) as employed in ref. [p2 |. Its 
straight-line input r = ru) = t — is with e > becomes strongly deformed after 
the transition to £ = £(r) = ((rw). 

In more detail, the left path could also inessentially be deformed, for the sake of 
simplicity, in such a way that e — e(t) — ► for \t\ — > oo. This asymptotically modified 
contour of r = tm avoids the upwards-running cut and comfortably coincides with 
the real line at both its asymptotic ends. It remains parametrized by t G (— oo, oo) 
and defines the right-hand side bent curve £ = Q — iZ via the pair of the real implicit 
equations 

sinhtcos£(t) = e z sinf2, 
cosh t sin e(t) = e z cosfl. 

Their explicit solution 



n 



(*) 



arctan 



tanh t 



tane(t) 



'-'(—oo) j ^(oo) 



-f +£(-00),f 



£ OO. 



'(*) 



\ In [sinh 2 1 + sin 2 e(t) 



shows that the new curve £ = = Q(t) — iZ(t) is arch-shaped and symmetric. It 
starts in a left imaginary minus infinity (at Z(_oo) ^ 1) an d ends in its right parallel 
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(with O^oo) = 7r/2 — e(oo) and ^> 1). Its top Z/ \ = In sine < at t = 0(o) = 
moves upwards in an inverse proportion to the decrease of the original shift e:(0) — > 0. 

III. PAIRS OF STATES WITH THE SAME N 

In a search for a "new" solvable model V^ D \^), let us start from the complexified 
Schrodinger equation ([2.1|) defined on a complex contour r = tm = t—ie(t). We shall 
apply the change of the complex contours ( |2.5| ) to the complexified Eckart potential 



V^r) = ^^-2i^, p>0. (3.1) 

smh r smhr 

Although our result will coincide (not surprisingly) with the above-mentioned shape- 
non-invariant potential (|1.2| ), the mapping itself will be shown to exhibit certain very 
unusual and unexpected features. 

In a preparatory step, the exact solvability of such a model may be most easily 
demonstrated via an auxiliary re-parametrization suggested in ref. ||22| , 

^ E \r) = {y- l)«(y + 1)> (i^) , y = y{r) = ^ = 1 - 2z 

with 

Av 2 = 2i(3 - E, Au 2 = -2i/3 - E. 
This leads to the new differential equation 

z(l-z) ip"(z) + [c - (a + b + l)z] ip'(z) - ab ip(z) = (3.2) 

which is completely solvable in terms of the Gauss hypergeometric series, 

tp{z) = Ci • 2 Fi{a, b; c; z) + C 2 ■ z l ~ c 2 F 1 (a + 1 - c, b + 1 - c; 2 - c; z) (3.3) 

where 

c = l + 2w, a + b = 2u + 2v + l, ab = (u + v)(u + v + 1) + A(l - A). 
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All the freedom of parameters becomes fixed by the boundary conditions ip(±oc) = 
which give b = —N and a = b ± {2 A — 1). Jacobi polynomials P^ a, ^{x) |23|] enter 
the elementary formula for the wave functions, 

^( £ )( r ) = ( _J_\ U+V e (v-u)r . ? W^) (cothr)i (3-4) 
V smh r / 

These solutions remain normalizable if and only if 

a = 2A- N -1, u + v = A-N-l>0, u - v = -i %■ — -. (3.5) 

A-N-l 

The model possesses N max < A — 1 bound states with the energies 

E (E) = _ {A _ N _ 1} 2 + ^_^_ i)2 1 N = 0, 1, ... , N max . (3.6) 

We are prepared to transform this solution into its partner of eq. ( |2.3|) . It suffices 
to change variables via eq. ( |2.5|) . This replaces the Eckart problem (|2.1| ) by the new 
Schrodinger equation fl2.3p . The implicit definition ( ^.4[ ) of the new potential \^ (D) (0 
acquires a more explicit form 



V (D) {r) _ E (D ) = ^_ 1 

4 cosh r cosh r 



N(N + 1) + (2N + l)S N + 1 + ^ 

-tuusii / UUS11 / L °N 

oi, / cosh r 



+ 



+ f^-4l+2^^^, 6 N = A-N-1>0. (3.7) 



In the domain of the large \r\ 3> 1 this formula is dominated by the last two terms. 
Only the very last one depends on the sign of Re r so that the coupling j3 must be 
independent of iV (we do not wish to have a state-dependent potential). We then 
determine (i.e., strictly speaking, remove the shift-ambiguity of) the energy E^ D > 
by the convenient requirement that V^ D \±oo) vanishes at (5 — 0. The other two 
asymptotically smaller components of V( D >(£) are of the first and second order in 
1/ cosh 2 r(£). Both of them must be also independent of iV of course. In the first 
order this gives the strict rule 

N 2 + N + 1 + (2N + 1) 5 + (3 2 /5 2 = constant (= C) . (3.8) 
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In the second order, the coefficient is equal to 3/4 and the condition remains trivial. 
We just confirmed that the replacement of r by r(£) transforms the Eckart potential 
into its exactly solvable descendant (|1.2|) , indeed. 
During the re-construction of the potential V^(^), energies 

4f ] and wave func- 
tions V'lv (£) the auxiliary, A^— dependent value of 5 = 5n = A — N — 1 > is to be 
determined as a root of the cubic equation (|3.8|). In Hermitian setting, the correct 



account of the physical boundary conditions makes this root unique [16|]. In the 
generalized, VT symmetric setting, the exceptional boundary condition in the origin 
becomes redundant. This is the reason why we have chosen our particular example. 
A priori, one may expect that the choice of the root 5n could be ambiguous. 

In the light of our present construction, the latter expectation proves fulfilled. At 
the sufficiently large values of C > C min > there exist three real roots Sn- Only 
one of them (viz., the negative one) can be eliminated as violating the asymptotic 
physical boundary conditions (i.e., the normalizability of the wave function). In 
contrast to the Hermitian case, two of the roots 5 = 8^ = 5^ (j3,C) > of our 
cubic eq. ( |3.8| ) remain equally acceptable. At any number of nodal zeros N, each of 
them defines a separate energy level, 

E = EffM =(s { N ±) + N + 1 -) 2 + 3 --C. (3.9) 



27 4 

The related VT symmetric wave functions 



^ = #(0 = fr(f)]-#[»-(f)] 



are proportional to the same Jacobi polynomials as above, 



ij} — e 



-i8£ 



1 - e- 2i « 



1/2 



This is the core of our present message. In the \t\ ^> 1 asymptotic regions we have 
£ ~ — iZ ± 7r/2 with Z ^> 1. This re-confirms that both our series of wave functions 
ip^ are asymptotically vanishing as exp (-6&Z). For both the roots 5 = 5^ > 0, 
they are safely normalizable. 
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IV. DISCUSSION 



The use of the commutativity [H, VT] = complies with our intuitive expec- 
tations in the numerical setting |24] and in perturbation theory |25 ], in the WKB 
approximation |26|] and in the supersymmetric context |27| as well as in the phe- 
nomenologically oriented field-theoretical studies I n contrast, our present results 
form a paradox since the change of variables mediates an utterly unusual one-to-two 
correspondence between the two complex potentials. For the force V^ D \^) this gives 
the two parallel series of bound states which must be distinguished by an additional, 
parity-type quantum number q — ±1. The nodal count iV itself does not suffice 



to characterize the energy levels This is a puzzling situation since we see 

no obvious reason for the introduction of q. One has to assume the existence of 
some unknown, hidden symmetry in our problem (|2.3| ), but it is still necessary to 
accept the fact that this symmetry has to break down during an apparently innocent 
change ( P~5| ) of the coordinates. 

We can summarize that our above conclusions are to be added to the list of 
the "oddities" which emerge in VT symmetric quantum mechanics due to it weaker 
boundary conditions. Let us just remind the reader that this is not in fact an 
isolated paradox. Our traditional intuition has already had really hard times with 
the unavoided level crossings in ref. [10], with the decrease of the energy with iV in 
ref. |20|| , with the high excitations caused by a weak potential in refs. [29] etc. Thus, 
just another item is provided by our present example. 

We may formulate our tentative conclusion that the VT symmetric deformations 
of the integration paths can in fact destroy the (a priori, plausible) similarity between 
the complex and ordinary parity. One of the most important implications is that 
any future appropriate generalization of the Hermitian Sturm-Liouville oscillation 
theorems will be necessarily not entirely trivial. Also the closely related concept of 
the completeness of states |30| must be dealt with an exceptional care in any future 
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development of the VT symmetric quantum mechanics. We might point out that 
the similar words of warning have been recently issued also on the purely numerical 
basis H3T . 
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